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We present a new method to compute quantum energies in presence of a back- 
ground field. The method is based on the string-inspired worldline approach to 
quantum field theory and its numerical realization with Monte-Carlo techniques. 
Our procedure is applied to the study of the Casimir force between rigid bodies 
induced by a fluctuating real scalar field. We test our method with the clas- 
sic parallel-plate configuration and study curvature effects quantitatively for the 
sphere-plate configuration. The numerical estimates are compared with the "prox- 
imity force approximation". Sizable curvature effects are found for a distance-to- 
curvature-radius ratio of a/R > 0.02. 



1 Introduction 

This work is based on the formulation of the Casimir effect in the framework 
of a renormalizable quantum field theory. The Casimir energy is accessed 
via the quantum effective action expressed in the string-inspired worldline 
approach to quantum field theory ^ Calculations are performed using a re- 
cently developed numerical method to compute quantum energies in presence 
of a background field El, which has already proved to be useful in the study of 
fermion-induced effects of magnetic interfaces^. In the next section, we de- 
rive the expression of the effective action on the worldline, define the Casimir 
energy and describe briefly the renormalization program. Sectionals devoted 
to the loop cloud method, our numerical implementation of worldline calcu- 
lations. Finally, we focus in section 0] on the study of Casimir forces between 
rigid bodies for the classic parallel-plate configuration, which serves as test for 
the method, as well as for the experimentally relevant sphere-plate configura- 
tion. Curvature effects are investigated and a comparison with the proximity 
force approximation is performed. 



2 Worldline formalism 



The embedding of the Casimir effect in quantum field theory, which is used in 
this approach and other contributions to this conference'^', is an alternative 
to the standard approach of imposing ab initio boundary conditions on the 
fluctuating field. By contrast, the quantum field theoretic approach permits a 
realistic modeling of the interactions with the physical boundary and a careful 
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study of the Casimir energy in the limiting case where boundary conditions 
are satisfied ^1 

We focus in this work on the Casimir effect induced by a fluctuating real 
scalar field. The starting point is given by the field theoretic lagrangian 

£ = \d^d^ + ^m 2 ^ 2 + \v{x)4, 2 . (1) 

Here we work in D = d + 1 Euclidean spacetime dimensions, i.e. d space 
dimensions. The interplay between the fluctuating field and the boundary 
matter is taken into account by the interaction term V(x)(j) 2 where the po- 
tential V(x) models the physical properties of the Casimir interface. The 
complete unrenormalized quantum effective action for V is 

m . ifc la - a2 \f + y. (2, 
1 1 2 -~d 2 + m 2 w 

The worldline representation of T[V] is obtained ® by (i) introducing the 
propertime representation of the functional logarithm with UV regularization 
(e.g. a cutoff A at the lower bound of the propertime integral), (ii) performing 
the trace in configuration space Tr[&\ — > J d D x(x\0\x) and (iii) interpreting 
the matrix element as a quantum mechanical transition amplitude expressed 
by a Feynman path integral. The final expression reads 

m = -— e-^J d°x CM [(W v [T,x(r)]) x - 1] , (3) 

where the expectation value is defined by 

(W v [T,x(r)]} x := X -^E1 . (4 ) 

/ Vx{t) e~Jo dT ~ 

x(0)=x(T) 

This expression is the central piece of the worldline approach to quantum 
field theory. From Eq. Q , we see that the computation of the effective action 
corresponding to the theory 0) with potential V is based on the average of the 
holonomy factor Wv[T, x(t)] over all worldlines x(t) which are parametrized 
by the propertime r £ [0, T] and closed, x(0) = x(T), due to the tracing 
operation in Eq.®. In Eq.@, all worldline loops have a common center of 
mass xcMj which implies J Q drx M (r) = 0. 

In this work, we concentrate exclusively on static Casimir configurations, 
i.e. the modeling potential 

V(x) = V(x) 

is time independent. In this case, the propertime integrand itself does 
not depend on time and the time integration can be carry out trivially, 
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J dxQ M — L Xo , where L Xo denotes the "volume" in time direction. We define 
the (unrenormalized) Casimir energy as 

£ V =T[V]/L X0 . (5) 

The effective action as given in Eq.Q is an alternative formulation of 
the quantum field theory defined from Q). For this reason, the divergence 
structure of the unrenormalized effective action on the worldline can be ex- 
pressed in the language of conventional Feynman diagrams. This is done by 
expanding the propertime integrand in Eq. for small values of T: 

J d D x(W v [T,x(T)} - l) x = -T J d D xV(x) + ^- J d D x V{xf + C(T 3 ), 

(6) 

which should be read together with the propertime factor \/T D / 2+1 of Eq.lJSJ. 
For D > 2, the first term is divergent and can always be eliminated by using 
the "no-tadpole" renormalization scheme. The second term is divergent for 
D > 4 and contributes to the one-loop diagram containing two insertions 
of the potential. Standard renormalization provides us with a counter term 
~ J dxV 2 subject to a physically chosen renormalization condition such that 
the divergence arising from the T 2 is canceled and which fixes the physical 
value of the renormalized operator ~ V 2 . 

This procedure permits to handle these divergences which have a clear 
field theoretic origin. Moreover further divergences may arise from the mod- 
eling of the Casimir configuration, i.e. from the potential V itself. These 
divergences occur when the potential is tuned to take some limiting form 
V — > V cr correspondingly to some convenient assumptions (ideal cases), 
r[V — > V ct ] — * oo, and may not be removed in a physically meaningful way. 
The relevant question then is as to whether physical observables are affected 
by these divergences or not. If not, the problem can possibly be bypassed. 



3 Loop Cloud Method 

As it is clear from Eq. © , the main task to be performed is the average (@J over 
the ensemble of all closed worldlines x(t) centered upon a given point xcm- 
The idea of the loop cloud method is to perform this expectation value by 
generating the loop ensemble in a numerical way^l. It is evident that we can 
generate neither the whole loop ensemble, which is infinite, nor an arbitrarily 
big amount of loops, for clear reasons of computational time. We will adopt 
a strategy which is widely used in statistical physics, namely the selection 
among the ensemble of all loops of the relevant configurations contributing 
to the loop average. These are generated according to the Gaufiian loop 

distribution pi 00 p[T, x(t)] ~ e Jo 4 . Generating an ensemble of N\ p such 
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loops {x(r)|i, i = 1, . . . , -/Vip}, the expectation value is approximated by 



(W v [T,x(T)}} x ~—J2Wv[T,x(T)\ i }. (7) 

The numerial burden is moreover dramatically reduced if we rescale the loops 
x M (r) in such a way that the distribution p\ oop becomes T-independent. Let 
us consider the following rcscaling transformations: 

r G [0, T] t := I G [0, 1], z(r) - y(t) := 

The propertime integral appearing in the loop distribution becomes 

dTX 2 { T )-> f dty 2 {t) 
Jo 

where the T-dependence has indeed disappeared. The numerical task is then 
reduced to the all-at-once generation ^ of an ensemble of y-loops, which are 
called unit loops. The evaluation of the expectation value 10 for a given 
value of T consists then simply in the rescaling of all members of the unit 
loop ensemble x(t) = \/Tv(t/T) before the sum (J7J) is performed. It is also 
obvious that the concept of continuous loops has to be abandoned. In the 
simulations a loop is defined by a collection of N pp i points, corresponding to 
the discretization of the propertime interval [0, T]: 

X ( T ) -* {^0 = x(t q ),x 1 = x(n), ... ,£at pp1 „i = x(rjv ppl _i),XAr ppl = x } . 

Let us point out that this discretization is not to be mistaken for a spacetime 
discretization. Every point Xi can move freely in continuous spacetime. 

The loop cloud method is developed independently of the background 
potential. This implies in particular that the Casimir effect can be studied 
numerically for arbitrary geometries and independently of the degree of sym- 
metry of the physical boundaries. This permits us to go beyond the geometry 
restrictions required by the 'proximity force approximation' (PFA)^™ which 
is the standard analytical approach to non trivial geometric configurations. 
Beyond this, the method does not require perfect geometries and perfect sur- 
faces but can, for instance, facilitate the study of corrugated surfaces^. 



4 Casimir forces between rigid bodies 

In this work, we model the physical boundaries by a potential of the form 

V(x) = A/(x), 

where /(x) describes the geometry of the Casimir configuration and A > 
is the strength of the field-matter interaction. It can roughly be viewed as 
a plasma frequency of the boundary matter: for fluctuations with frequency 
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lu 3> A, the Casimir boundaries become transparent. Two types of limits can 
be taken regarding this parametrization: 

/(x) — > / da S d (x — x CT ) (sharp potential limit), (8) 
Jt, 

A —* oo (strong coupling limit). (9) 

Here £ represents the geometry of the Casimir configuration and denotes a 
d — 1 dimensional surface, generally disconnected (e.g. two separated plates, 
£ = Si + S2). The first limit idealizes the plate to be infinitely thin. The sec- 
ond limit imposes the Dirichlet boundary condition, implying that all modes 
have to vanish on S. Together they correspond to the case of ab initio bound- 
ary conditions of the standard approach. These limits are unfortunately 
spoiled by the occurence of a divergence in the Casimir energy which can- 
not be removed by the standard field theoretic renormalization procedure. 
However, this problem can be bypassed if we consider instead the Casimir 
force acting on the rigid bodies. The force is defined by 

-Fcas(a) = --^-E Cas (a) 
oa 

where a denotes the distance between the boundaries. This definition gives 
us a certain freedom of the choice of the Casimir interaction energy i?cas («) ■ 
For a Casimir configuration of two rigid bodies whose surfaces are denoted by 
Si and 5*2, we define the interaction energy as 

E Cas (a) = £v 1+2 (a) - £v x - £v 2 (10) 

where £y 1+2 , £vi and £v 2 are the energies (jSJ) with the potentials V1+2, Vi and 
V2 given in the limiting form JSJ) where £ = SiU S2, £ = Si and £ = S2, 
respectively. The first term depends on the distance a and contributes to the 
Casimir force. Since the remaining terms do not depend on a, the way there 
are chosen has no influence on the Casimir force. In the present case, they 
are chosen in such a way that -Ecas (a) is rendered finite by substracting each 
Casimir energy of the single bodies. 

4-1 A benchmark test: the parallel plates revisited 

Casimir forces can be analytically computed only for a small number of rigid- 
body geometries among which there is the parallel-plate configuration. Com- 
paring our numerical estimates with the analytically known result ^| the 
agreement is very satisfying in all regimes of scalar mass m, coupling A and 
distance a , see Fig^ right panel. Moreover, we have found that our algo- 
rithm is scalable: if higher precision is required, only the parameters of the 
loop ensemble, points per loop N pp \ which controls the systematic error in- 
duced by the loop discretization and number of loops N\ p which controls the 
approximation of the loop average l@J by the sum 0, have to be adjusted, 
as seen in FigQ] left panel. 
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Figure 1. Numerical estimate of the interaction Casimir energy per unit area of the parallel- 
plate configuration for various loop ensembles as a function of the number of points per 
loop N pp i (left panel). Numerical estimates vs. analytical results as a function of the 
dimensionless quantity ma, A = 100m (right panel). 



4-2 Beyond the PFA: the sphere-plate configuration 

We study the experimentally relevant sphere-plate configuration and confront 
our numerical estimates with the results provided by the proximity force ap- 
proximation (PFA). Due to the assumptions on which the PFA is based, i.e. 
dropping of non-local curvature effects, this approach is expected to give rea- 
sonable results only if the typical curvature radii of the surfaces elements 
is large compared to the element distance. Beside its physical interest, the 
sphere-plate configuration with its relatively simple geometry provides us with 
a beautiful illustration of the loop cloud machinery. As illustrated in Fig|2 
the main feature of the sphere-plate physics, i.e. curvature effects, can be 
understood qualitatively simply by thinking in terms of loops. On the one 
hand, the energy density at a point xq is obtained by performing averages 
over loop clouds centered on xq. On the other hand, it can be shown that the 
energy density is proportional to the number of loops piercing both physical 
boundaries. With these two statements only, Fig.(0) can be qualitatively un- 
derstood. The effect of the increasing curvature of one of the boundaries is 
the decrease of contributing loops near the plate, implying the decrease of the 
Casimir energy in this region. 

Let us finally consider the complete interaction Casimir energy for the 
sphere-plate configuration as a function of the sphere-plate distance a (we 
express all dimensionful quantities as a function of the sphere radius R). In 
Fig. |21 we plot our numerical results in the range a/R ~ 0(0.001 . . . 10). Since 
the energy varies over a wide range of scales, already small loop ensembles 
with rather large errors suffice for a satisfactory estimate (the error bars of 
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Figure 2. Interaction Casimir energy density along the symmetry axis (x axis) for the 
sphere-plate configuration in comparison to the parallel-plate case. Close to the sphere, the 
worldline loops do not "see" the curvature; but at larger distances, curvature effects enter 
the energy density. For illustration, the sphere-plate geometry is also sketched (thin black 
lines) . 



an ensemble of 1500 loops with 4000 ppl cannot be resolved in Fig. EJ) • Let 
us compare our numerical estimate with the proximity force approximation 
with radial distance measurements: using the plate surface as the integration 
domain for the PFA, we obtain the solid line in Fig. (PFA, plate-based), 
corresponding to a "no-curvature" approximation. As expected, the PFA ap- 
proximation agrees with our numerical result for small distances (large sphere 
radius). Sizable deviations from the PFA approximation of the order of a 
few percent occur for a/R ~ 0.02 and larger. Here, the curvature-neglecting 
approximations are clearly no longer valid. This can be read off from the 
right panel, where the resulting interaction energies are normalized to the 
numerical result. 

5 Conclusion 

In this work, we have developed a new numerical method to compute quan- 
tum energies in presence of a background field. The algorithm is constructed 
independently of the background, which opens a wide range of applicability 
of the procedure. In the context of the Casimir effect, arbitrary geometric 
configurations can be considered. We focused in this work on the Casimir 
forces induced on rigid bodies by the fluctuations of a real scalar field. The 
method was tested by comparing our numerics with the analytical results in 
the classic parallel-plates configuration and confronted with the results pro- 
vided by the proximity force approximation in the sphere-plate configuration. 
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Figure 3. Left panel Logarithmic plot of the interaction Casimir energy for the sphere- 
plate configuration. For small separations/large spheres, a/R < 0.02, the proximity force 
approximation (PFA) approximates the numerical estimate well; but for larger a/R, cur- 
vature effects are not properly taken into account. The PFA becomes ambiguous for larger 
a/R, owing to possible different choices of the integration domain in the PFA. A geomet- 
ric mean (dotted-dashed line) of the limiting cases shows reasonable agreement with the 
numerical result. Right panel Interaction Casimir energies normalized to the numerical 
result . For a/R > 0.02, the fluctuation-induced curvature effects occur at the percent level. 



We have been able to study the usually neglected nonlocal curvature effects 
which become sizable for a distance-to-curvature-radius ratio of a/R > 0.02. 
Some generalizations to more realistic configurations are straightforward: fi- 
nite temperature effects, surface roughness and corrugation. Let us finally 
note that the implementation of worldline numerics for the Casimir effect due 
to a fluctuating electromagnetic field has still to be worked out. Here the new 
ingredient is to find an efficient field theoretic formulation of the interaction 
of the electromagnetic field with the Casimir boundaries. 

Acknowledgments 

L.M. is grateful to the organizers of the QFEXT03 meeting in Norman for the 
stimulating atmosphere of the conference and acknowledges financial support 
by the Deutsche Forschungsgemeinschaft under contract GRK683. H.G. is 
supported by the Deutsche Forschungsgemeinschaft under contract Gi 328/1- 
2. 



References 

1. H. Gies, K. Langfeld and L. Moyaerts, JHEP 0306 (2003) 018. 

2. H. Gies and K. Langfeld, Nucl. Phys. B 613 (2001) 353, 

H. Gies and K. Langfeld, Int. J. Mod. Phys. A 17 (2002) 966. 



8 



3. K. Langfeld, L. Moyaerts and H. Gies, Nucl. Phys. B 646 (2002) 158. 

4. N. Graham, R. L. Jaffe, V. Khemani, M. Quandt, M. Scandurra and 
H. Weigel, Phys. Lett. B 572 (2003) 196. 

5. R. L. Jaffe, |arXiv:hep-th/030foT^ 

6. For a review see C. Schubert, Phys. Rept. 355 (2001) 73. 

7. B.V. Derjaguin, I.I. Abrikosova, E.M. Lifshitz, Q.Rev. 10 (1956) 295. 

8. J. Blocki, J. Randrup, W.J. Swiatecki, C.F. Tsang, Ann. Phys. (N.Y.) 
105 (1977) 427. 

9. T. Emig, Europhys. Lett. 62 (2003) 466. 

10. M. Bordag, D. Hennig and D. Robaschik, J. Phys. A 25 (1992) 4483. 



9 



